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^ The circumgalactic medium (CGM) is fed by galaxy outflows and accretion 

^ of intergalactic gas, but its mass, heavy element enrichment, and relation to 

galaxy properties are poorly constrained by observations. In a survey of the 
>^ outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the 

^ Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos 

^ of ionized oxygen surrounding star-forming galaxies, but we find much less 

ionized oxygen around galaxies with little or no star formation. This ionized 
CGM contains a substantial mass of heavy elements and gas, perhaps far ex- 
ceeding the reservoirs of gas in the galaxies themselves. It is a basic component 
of nearly all star-forming galaxies that is removed or transformed during the 
quenching of star formation and the transition to passive evolution. 

Galaxies grow by accreting gas from the intergalactic medium (IGM) and converting it to 
stars. Stellar winds and explosions release gas enriched with heavy elements (or metals, i), 
some of which is ejected in galactic-scale outflows (2). The circumgalactic medium (CGM), 
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loosely defined as gas surrounding galaxies within their own halos of dark matter (out to 100- 
300 kiloparsec), lies at the nexus of accretion and outflows, but the structure of the CGM and 
its relation to galaxy properties are still uncertain. Galactic outflows are observed at both low 
(2-4) and high (5-7) redshift, but it is unclear how far they propagate, what level of heavy- 
element enrichment they possess, and whether the gas escapes the halo or eventually returns 
to fuel later star formation. Models of galaxy evolution require significant outflows to explain 
observed galaxy masses and chemical abundances and to account for metals observed in the 
more diffuse IGM (8, 9). The CGM may also reflect the theoretically -predicted transition from 
filamentary streams of cold gas that feed low mass galaxies to hot, quasi-static envelopes that 
surround high mass galaxies (10, 11). Both outflow and accretion through the CGM may be 
intimately connected to the observed dichotomy between blue, star-forming, disk-dominated 
galaxies and red, passively evolving, elliptical galaxies with little or no star formation (72). 
However, the low density of the CGM makes it extremely difficult to probe directly, thus models 
of its structure and influences are typically constrained indirectly by its effects on the visible 
portions of galaxies, not usually by observations of the gas itself. 

We have undertaken a large program with the new Cosmic Origins Spectrograph (COS) 
aboard the Hubble Space Telescope to directly map the CGM using the technique of absorption- 
line spectroscopy, in which a diffuse gas is detected by its absorption of light from a background 
source. Our background sources are UV-bright quasi-stellar objects (QSOs), which are the 
luminous active nuclei of galaxies lying far behind the galaxies of interest. We focus on the 
ultraviolet 1032,1038 A doublet of O VI (0+^), the most accessible tracer of hot and/or highly 
ionized gas at redshift 2; < 0.5. O VI has been used to trace missing baryons in the IGM (14-17), 
the association of metals with galaxies (18-20), and coronal gas in the Milky Way halo (20). 

The high sensitivity of COS enables a QSO absorption-line survey of halos around galaxies 
with a predetermined set of properties. We have selected 42 sample galaxies (Tables SI and 
S2) that span redshifts Zgai = 0.10 - 0.36, stellar masses log(M*/MQ) = 9.5 - 11.5. The 
QSO sightlines probe projected radial distances (i.e. impact parameters) to the galaxies of 
i? = 14 — 155 kpc. We used the COS data to measure the O VI column densities (iVovi 
in cm~^), line profiles, and velocities with respect to the target galaxies (Fig. 1, 21). We 
measured precise redshifts, star formation rates (SFR in Mq yr~^), and metallicities for our 
sample galaxies using low-resolution spectroscopy from the Keck Observatory Low-Resolution 
Imaging Spectrograph (LRIS) and the Las Campanas Observatory Magellan Echellette (MagE) 
spectrograph (22, SOM text S2). 

Our systematic sampling of galaxy properties allows us to investigate the connection be- 
tween galaxies themselves and the CGM. The O VI detections extend to R = 150 kpc away 
from the targeted galaxies, but the whole sample shows no obvious trend with radius R (Fig. 
2). The strong clustering of detections within ± 200 km s~^ of the galaxy systemic velocities 
indicates a close physical and/or gravitational association. 

Circumgalactic medium gas as traced by O VI reflects the underlying bimodality of the 
general galaxy population (12, 23). We find a correlation of Nqyi with specific star formation 
rate sSFR (= SFR/M*; Fig 3). For the 30 galaxies with sSFR > 10~^^ yr-\ there are 27 
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detections with a typical column density logAi'ovi = 14.5 {24) and a high covering fraction 
/hit ~ 0.8 — 1 maintained all the way to i? = 150 kpc (Fig. 2). For the 12 galaxies in the passive 
subsample (sSFR < 10^^^ yr^^) there are only 4 detections with lower typical Nqyi than the 
star- forming subsample (25). Accounting for the upper limits in A^'ovi and sSFR, we can reject 
the null hypothesis that there is no correlation between A^ovi and sSFR at > 99.9% confidence 
for the whole sample and > 98% for each of the 50 kpc annuli shown in Figure 1 (SOM text 
S5). This effect remains even when we control for stellar mass: a Kolmogorov-Smirnov test 
over log Af^, > 10.5, where the star-forming and passive subsamples overlap, rejects at > 99% 
confidence the null hypothesis that they draw from the same parent distribution of Nq^i (Fig- 
S2). We therefore conclude that the basic dichotomy between star-forming ("blue-cloud") and 
passive ("red- sequence") galaxies is strongly reflected in their gaseous halos, and that the CGM 
out to at least 150 kpc either directly influences or is directly affected by star formation. 

O VI is a fragile ionization state, which never exceeds a fraction /ovi = 0.2 of the total 
oxygen for the physical conditions of halo gas and is frequently much less abundant (Fig. 4). 
Our observations imply a typical CGM oxygen mass for star-forming galaxies of 

Mo = 57ri?2(iVovi)mo/hit \^] = 1-2 x 10^ ( ^ ) Mq , (1) 

V /OVI / V Jovi / 

where we have taken a typical {Nqvi} = 10^^'^ cmr"^ and R = 150 kpc, and the hit rate cor- 
rection /hit computed separately in three 50 kpc annuli (Figs. 1 and 2). This mass of oxygen 
is strictly a lower limit because we have scaled to the maximum /ovi = 0.2 (Fig. 4). The 
corresponding total mass of circumgalactic gas is 

M,.. = 177 Mo = 2 X 10^ (^) Me , (2) 

where Z is the gas metallicity, and the solar oxygen abundance is ^^o/^^H = 5 x 10^^ (2(5). 

Even for the most conservative ionization correction (/ovi = 0.2), the OVI-traced CGM 
contains a mass of metals and gas that is substantial by comparison with other reservoirs of 
interstellar and circumgalactic gas. If our sample galaxies lie on the mean trend of gas frac- 
tion for low-z galaxies (27), they have interstellar medium (ISM) gas masses of Mjsm = 
(5 — 10) X IO^Mq and contain Mjgp^ = (2 — 10) x lO^Af© of oxygen, taking into account 
the observed correlation between galaxy stellar mass and ISM metallicity (Fig. 4, SOM text 
S5). The minimum CGM oxygen mass is thus 10-70% of the ISM oxygen (Figs. 4 and S4). 
The covering fractions and column densities we find for star-forming galaxies are insensitive 
to M^,, while the ISM metal masses decline steeply with M* according to the mass-metallicity 
relation. Thus the ratio of CGM metals to ISM metals appears to increase for lower mass 
galaxies (assuming constant fovi), perhaps indicating that metals more easily escape from their 
shallower gravitational potentials. The implied total mass of circumgalactic gas Mgas is more 
uncertain because it can strictly take on any metallicity; for a fiducial solar metallicity. Equa- 
tion (|2]) implies a total CGM mass comparable to Mism, and several times the total mass inferred 
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for Milky Way "high velocity clouds" (HVCs) (28, 29) or for low-ionization (Mg II) gas filling 
halos toR= 100 kpc {30). 

For the densities typically expected at radii R ~ 100 kpc, /ovi exceeds 0.1 only over a 
narrow temperature range lo^ "^"^-^ K, and it only exceeds 0.02 over lo^-^"^-^ K (Fig. 4). Either 
a large fraction of CGM gas lies in this finely tuned temperature range — a condition that is 
difficult to maintain because gas cooling rates peak at T fa 10^^ K — or the CGM oxygen and 
gas masses are much larger than the minimum values we have quoted above. Lower density, 
photoionized gas can achieve high /ovi ~ 0.1 over a wider temperature range, but at these low 
densities it is hard to produce a 10^^'^ cm~^ column density within the confines of a galactic 
halo, especially if the metallicity is low (Fig. S5). Thus /ovi = 0.02 and Z — O.IZq are 
plausible conditions for the O Vl-traced gas, but it is unlikely that both conditions hold simul- 
taneously. However, if either one holds the CGM detected here could represent an important 
contribution to the cosmic budgets of metals and baryons. In either case Mgas is compara- 
ble to the total ~ 3 x lO^^M© inside R = 300 kpc inferred from H I measurements at low 
redshift (19), and the ~ 4 x lO^^M© inferred for the CGM surrounding rapidly star-forming 
galaxies at 2; fa 2 — 3 (31). By generalizing our typical Mq to all star-forming galaxies with 
M* > 10^-^ Mq, we estimate that the halos of such galaxies contain 15% x (0.02//ovi) of the 
oxygen in the universe and 2% x (0.02//ovi) x (Zq/Z) of the baryons in the universe. 

The metals detected out to i? ~ 150 kpc must have been produced in galaxies, after which 
they were likely transported into the CGM in some form of outflow. However, these outflows 
need not be active at the time of observation; indeed the large masses imply long timescales. 
Because 1 Mq of star formation eventually returns 0.014 Mq of oxygen to the ISM (32), at least 
8.6 X IO^A/q of star formation is required to yield the detected oxygen mass. This is equivalent 
to ~ 3 X 10^ yr of star formation at the median SFR = 3 Mq yr"^ of our star-forming sample, 
in the unlikely event that all oxygen produced is expelled to the CGM, and longer in inverse 
proportion to the fraction of metals retained in the ISM. Thus the detected oxygen could be the 
cumulative effect of steady enrichment over the preceding several Gyr, the product of sporadic 
flows driven by rapid starbursts and an active nucleus (33), or the fossil remains of outflows 
from as early as z ^ 1.5 — 3 (7, 31). While the exact origin of the mass-metallicity relation of 
galaxies is not yet known, models that explain it in terms of preferential loss of metals imply 
that a substantial fraction of the metals produced by star formation must be ejected from the 
galaxy rather than retained in the ISM (27). The CGM detected here could be a major reservoir 
of this ejected material, with important consequences for models of galactic chemical evolution. 

The O VI we observe arises in bulk flows of gas over 100 — 400 km s~\ but the relative 
velocities are usually below halo escape speeds (Fig. 2), even when we take projection effects 
into account (Fig. SI). Thus much of the material driven into the halo by star formation could 
eventually be reacquired by the galaxy in "recycled winds," which may be an important source 
of fuel for ongoing star formation (34). It is unlikely that the detected gas is predominantly 
fresh material accreting from the IGM because models of "cold mode" accretion predict very 
low metallicity and low covering fractions /hit ~ 10 — 20% (35, 36), and "hot mode" accretion 
typically involves gas at temperatures T > 10^ K with undetectably low /ovi- 
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The passive galaxies in our sample once fonned stars, thus it follows that they would once 
have possessed hales of ionized, metal-enriched gas visible in O VI. The relative paucity of O VI 
around these galaxies implies that this material was transformed by processes that plausibly 
accompany the quenching of star formation (37), such as tidal stripping in group environments, 
re-accretion onto the galaxy in ionized form, or heating or cooling to a temperature at which 
O VI is too rare to detect. Our findings present a quantitative challenge for theoretical models 
of galaxy growth and feedback, which must explain both the ubiquitous presence of massive, 
metal-enriched ionized halos around star-forming galaxies and the fate of these metals after star 
formation ends. 
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Figure 1: An illustration of our sampling technique and data. (A) An SDSS composite 
image of the field around the QSO J1016+4706 with two targeted galaxies, labeled Gl and 
G2, which are both in the star-forming subsample. (B) The complete COS count-rate spectrum 
(counts s^^) versus observed wavelength. The panels below concentrate on the redshifted O VI 
1032,1038 A doublet for galaxies Gl (C) and G2 (D). The upper right panels illustrate the full 
sample by showing the locations of all sightlines in position angle and impact parameter R with 
respect to the targeted galaxies, for the star-forming (E) and passively evolving (F) subsamples. 
The circles mark R = 50, 100, and 150 kpc. 
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Figure 2: O VI association with galaxies. (A) O VI column density, A/'ovi in cm^^, vs. R 
in kiloparsec for the star-forming (blue) and passive (red) subsamples. Filled and open sym- 
bols mark O VI detections and 3cr upper limits, respectively. The detections in the star-forming 
galaxies maintain log A^'ovi ~ 14.5 to i? ~ 150 kpc, the outer limit of our survey. (B) Com- 
ponent centroid velocities with respect to galaxy systemic redshift for O VI detections, versus 
inferred dark-matter halo mass. The range bars mark the full range of O VI absorption for each 
system. The inset shows a histogram of the component velocities. The dashed lines mark the 
mass-dependent escape velocity at i? = 50, 100, and 150 kpc from outside to inside. 
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Figure 3: O VI correlation with galaxy properties. (A) O VI column density versus sSFR 
(= M^/SFR). Star- forming galaxies are divided from passively evolving galaxies by sSFR 
~ 10~^^ yr^^, our detection limit is sSFR ~ 5 x 10'^^ yr~^)- (B) The galaxy color-magnitude 
diagram (sSFR versus M J for SDSS-i-GALEX galaxies from (23). 
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Figure 4: CGM oxygen masses compared to galactic reservoirs. (A) The curves and the 
axis labels at right show the fraction of gas-phase oxygen in the O VI ionization state /ovi 
as a function of temperature, for three overdensities relative to the cosmic mean, p/p. All 
p/p> 1000 track the black curve on which coUisional ionization dominates, while for lower 
values photoionization by the extragalactic background can increase /ovi at low T. For gas 
that traces dark matter, p/p = 1000 is typical at i? ^ 100 kpc, while p/p = 50 — 100 for the 
outskirts of the halo. The pale green band shows the expected oxygen mass of the galaxies' 
ISM if they lie on the standard relation between Mjsm and and follow the mass-metallicity 
relation (MZR). The green dashed line shows the oxygen mass produced by 3 x 10^ Mq of star 
formation. The yellow band shows the expected oxygen mass for the extreme assumption that 
the typical host dark-matter halos (2 x 10^^ — IO^^Mq) have the universal baryon fraction and 
solar metallicity. (B) The CGM oxygen masses compared with the interstellar oxygen mass as 
a function of M^:. Points with range bars show the CGM oxygen mass Mq implied by Equation 
(1) for /ovi = 0.2, calculated separately for star-forming (blue) and passive (red) galaxies 
according to the hit rates in four bins of stellar mass. The purple curves show the calculated 
Mjgj^ for typical star forming galaxies in the SDSS, accounting for the mean MZR in the central 
curve and its uncertainties in the shaded region. The data points increase their mass in inverse 
proportion to /ovi- 
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SI Survey Design and Data Collection 

In recent years, the joint advances of the SDSS DR5 QSO catalog (40), the GALEX all-sky 

UV survey (41), and the installation of the Cosmic Origins Spectrograph (42) on-board HST 
have increased the numbers of UV-bright QSOs that can be observed with HST by an order of 
magnitude or more. We can now select a sample of QSOs whose sightlines pierce the halos 
of foreground galaxies with a predetermined range of properties, thus establishing a uniform 
sampling of galaxy stellar mass (M*), star-formation rate (SFR), and impact parameter (R). 

We began our sample selection by choosing QSOs atz<l that were bright enough to obtain 
COS spectra with S/N ~ 8 — 10 in five HST orbits or less. We require Zga_i > 0.11 to place the 
O VI doublet within the A > 1140 A wavelength range of COS. From the hst of approximately 
1000 such candidates that met all these criteria, we constructed an HST program of 39 QSOs 
and 134 orbits, with the galaxies sampling the range log M^/Mq = 10 — 11 (with some scatter 
out of this original range caused by errors in photometric redshifts) and i? = — 150 kpc as 
uniformly as possible. 

We obtained COS FUV spectroscopy of this sample of QSOs between March and Octo- 
ber 2010. The spectra were obtained at a resolving power of AA/A ~ 18,000—20,000 over 
1140 - 1750 A using the G130M and G160M gratings of COS. After standard processing by 
the CALCOS pipeline (v2.12), the individual exposures were co-added in counts space to pro- 
vide a final, reduced ID spectrum trace (using a method detailed by (43, 44)). 

Out of the 39 program targets, six COS spectra either did not provide coverage of the O 
VI doublet with data of acceptable quality, or they require special processing that is not yet 
implemented within CALCOS. We fully analyzed a total of 33 program datasets for this study. 
Our final sample consisted of 42 galaxies in the foreground of these QSOs: 30 star-forming and 
12 passively evolving galaxies. Eleven of the galaxies were not part of the original selection, 
although they meet the original criteria. These "bonus" galaxies were confirmed at different 
redshifts along the line of sight by optical spectroscopy, and we have checked that these galaxies 
do not bias the overall sample in any of the galaxy properties. A full presentation of the galaxy 
survey technique, the optical spectroscopy, and the associated errors is presented by (22). We 
exclude from the present sample those galaxies that did not provide good coverage of O VI in the 
COS data; the 42 galaxies in the present sample are listed in Tables SI and S2 for star-forming 
and passive galaxies, respectively. We label our galaxies with two numbers: the position angle 
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from the QSO, N through E, and the angular separation in arcseconds. The QSO and this galaxy 
designation can be used to cross-reference to the tables of the larger database in (22). 

To measure the precise redshifts and star formation rates of the galaxies, we obtained spectra 
of the galaxies with the Low-Resolution Imaging Spectrometer (LRIS) on the Keck telescope, 
and the Magellan Echellette (MagE) spectrograph on the Magellan Clay Telescope. With LRIS, 
we employed a 1" slit and the D560 dichroic, using the 600/4000 1/mm grism (blue side) and 
600/7500 1/mm grating (red side). This setup resulted in spectral coverage between 3000 and 
5500 A (blue side), and 5600 to 8200 A (red side). On the blue side, binning the data 2x2 
resulted in a dispersion of 1.2 A per pixel and a FWHM resolution of '--'280 km/s. On the red 
side, the data were binned 1x2, resulting in a dispersion of 2.3 A per pixel and a FWHM 
resolution of ~200 km/s. Exposure times were typically 800s in the blue and 2 x 360 s in the 
red, which resulted in signal-to-noise ratios of at least 3 per pixel for strong nebular emission 
lines in the galaxy spectra. Data reduction and calibration were carried out using the LowRedux 
IDL software package (45), which includes flat fielding to correct for pixel-to-pixel response 
variations and larger scale illumination variations, wavelength calibration, sky subtraction, and 
an initial flux caUbration using the spectrophotometric standard star G191B2B. 
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Table SI: O VI Sample Galaxies (Star-forming) 



QSO 


Galaxy 


iVovi cm-2 


J0042-1037 


358_9 


14.78 ± 0.08 


J040 1-0540 


67^4 


14.57 ±0.03 


J0820+2334 


260_17 


14.43 ±0.12 


J0910+1014 


34_46 


14.70 ±0.10 


J09 14+2823 


41^7 


14.69 ± 0.04 


J0925+4004 


193_25 


< 13.94 


J0943+0531 


106 34 


< 13.56 


J0943+0531 


227_19 


14.44 ±0.23 


J1009+0713 


170 9 


15.07 ±0.02 


J1009+0713 


204_17 


15.00 ±0.03 


J1016+4706 


274_6 


14.86 ± 0.02 


J1016+4706 


359_16 


14.70 ±0.03 


Jl 112+3539 


236_14 


14.70 ± 0.04 


Jl 133+0327 


164_21 


14.43 ±0.08 


J1233+4758 


94_38 


14.38 ± 0.06 


J1233-0031 


168_7 


14.68 ±0.03 


J1241+5721 


199_6 


14.81 ±0.02 


J1241+5721 


208 _27 


14.76 ± 0.02 


J1245+3356 


236_36 


14.34 ± 0.05 


J1330+2813 


289_28 


14.38 ±0.04 


J1342-0053 


157_10 


14.48 ±0.03 


J 14 19+4207 


132_30 


14.59 ±0.05 


J1435+3604 


68_12 


14.73 ± 0.04 


J1435+3604 


126_21 


14.57 ±0.05 


J1437+5045 


24_13 


14.51 ±0.12 


J 1437+5045 


317_38 


14.57 ±0.18 


J 1445+3428 


232_33 


14.29 ± 0.06 


J1550+4001 


97_33 


< 13.68 


J1555+3628 


88_11 


14.56 ±0.05 


J1619+3342 


113_40 


14.27 ±0.05 
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Table S2: O VI Sample Galaxies (Passive) 



QSO 


Galaxy 


NoYi 


J0226+0015 


268_22 


< 13.53 


J0910+1014 


35_14 


14.21 ±0.18 


J0928+6025 


110_35 


< 14.19 


J0935+0204 


15_28 


< 13.96 


J0943+0531 


216_61 


< 13.80 


J0950+4831 


177_27 


14.28 ± 0.04 


Jl 157-0022 


230 7 


< 13.56 


J1220+3853 


225 _38 


14.37 ±0.07 


J1342-0053 


77 10 


< 13.38 


J1550+4001 


197_23 


13.73 ±0.13 


J1617+0638 


253 39 


< 14.17 


J2257+1340 


270_40 


< 14.14 



The MagE spectrograph provides moderate resolution spectral coverage between 3200 A and 
10000 A. These data were acquired with a 0.7" slit and binned 1x1, giving typical dispersions 
of 0.3 A per pixel at 4500 A, and 0.5 A per pixel at 7800 A. The average emission-line FWHM 
is 55 km s~^. Exposure times varied between 600 and 1200 seconds, depending on the target 
galaxy apparent magnitude. The MagE spectra were reduced using the MASE pipeline (46). 
One-dimensional spectra are optimally extracted from the 2D reduced images. Spectrophoto- 
metric standard stars taken at a variety of airmasses were used to flux-calibrate the data. 

We make three corrections to the flux scale of our optical spectra: 1. A correction for slit 
losses that brings the spectra to an absolute flux scale. 2. A correction for Galactic extinction, 
assuming the E{B — V) value from (47) and a Galactic extinction law (48). 3. A correction for 
interstellar reddening to all line measurements from the observed H7 to H/3 and Ha to H/3 ra- 
tios. To bring the LRISb spectra to an absolute flux scale, we convolve it with the SDSS ^f-band 
filter response curve and scale the flux of the blue spectrum to match the reported SDSS-DR7 
petrosian magnitude. Similarly, we match the i-band magnitude for our LRISr spectra. The me- 
dian values of these two flux factors are 1.94 (blue) and 1.77 (red), corresponding to median slit 
light losses of 48% and 43%. These corrections, which provide accurate (at least as accurate as 
SDSS photometry) absolute fluxing, account for both slit losses and red/blue side matching in 
a self -consistent global fit. Furthermore, when performing B aimer corrections, we use as many 
Balmer lines as are detected in the spectrum and find generally good agreement after we have 
corrected the red and blue sides for slit losses and differences in absolute flux scales on the red 
and blue sides. However, sometimes there are only two Balmer lines detected, and occasion- 
ally, only one. For Balmer emission-line measurements, we account for contamination from 
underlying stellar absorption by fitting the continuum in the trough of detectable absorption. 
The overall effect on the line flux of the Balmer absorption ranges from 10% to 60% for the H^ 
emission line (when absorption is apparent). As for the LRIS data, we account for slit-losses 
by scaling the spectra to match SDSS r-band photometry (MagE contains no dichroic, so only 
a single band is used). 
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S2 Galaxy Properties from Spectroscopy 



Precise and accurate systemic redshifts were obtained for galaxies using a modified version 
of the SDSS IDL code "zfind", which works by fitting smoothed template SDSS eigenspectra 
to the galaxy emission-line spectra on both the red and blue sides. To account for systematic 
uncertainties in the absolute wavelength calibration and instrument flexure during the LRIS 
exposures, we adopt a conservative 30 km s~^ error on the galaxy redshifts. 

To obtain an estimate of the stellar mass of each galaxy, we used version 4_2 of the kcorrect 
IDL package {49), the SDSS DR7 galactic reddening corrected, asinh ugriz magnitudes, and 
the zfind spectroscopic redshifts. The stellar masses output by the routine sdssJccorrect are 
computed using a Chabrier IMF with the 1994 Padova isochrones. We raise the masses by a 
factor of 1.65 (0.22 dex) to compare them with quantities determined using a Salpeter IMF 
(i.e. SFRs, described below) and assign errors of 0.2 dex on propagated from the SDSS 
photometry and the systematic errors in the adopted mass conversion {50, 51). The masses 
we adopt here have been corrected assuming the 5 -year WMAP cosmology with a Hubble 
parameter h = 0.72 (52). 

We then calculate a current SFR using the Balmer emission lines Ha and H/3. For the former, 
we use the calibration of {53) where SFR [Mq yr'^] = 7.9 x 10"^^ Lna [ergs s"^]. This SFR 
calibration is derived from stellar population synthesis models that assume a Salpeter IMF and 
solar metallicity. When Ha is not observed in a spectrum, we use the same SFR calibration (55) 
divided by a factor of 2.86 for the H/3 emission line. This factor of 2.86 is the intrinsic ratio of 
Ha/H/3 at an effective temperature of 10,000 K and electron density of 100 cm^"^ for Case B 
recombination {54). Thus, the SFRs derived from Ha and H/3 are identical when we calculate 
a dust correction that gives Ha/HP = 2.86. The SFRs used in comparisons in this paper adopt 
the measurement from the best available line, or both as the case may be. 

When a galaxy's spectrum contains no emission lines, we measure an upper limit to the SFR 
by measuring the boxcar noise at the positions of H/3 and Ha. We use 3a line flux limits as our 
SFR upper limits in these cases, approximately 1/3 of our sample. We adopt conservative 3a 
limits to the SFRs since in these galaxies we are unable to make a correction for dust. 

S3 O VI Line Identification and Classification 

We begin by shifting the reduced, co-added COS spectra to the rest-frame of the target 
galaxy, setting v = Oto its measured systemic redshift. We then generate apparent optical depth 
plots (55) for the two lines of the O VI doublet to simplify their identification and confirmation, 
and to easily identify contaminating absorption. We then inspect the data to look for O VI 
absorption within ±600 km s^^ of the systemic redshift. Most of the O VI absorbers are readily 
apparent and easily confirmed in both lines of the doublet, but the redshifts of some targeted 
systems place one line of the doublet where it cannot be detected, either off the 1 140 A short 
end of the COS detector (zovi < 0.1), or underneath Galactic Lya (2;ovi ~ 0.17) or geocoronal 
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O I airglow emission (^ovi ~ 0.25). 

Contaminating absorption from other intervening absorption-line systems, or from the fore- 
ground ISM of the Milky Way, is always a risk. We attempt to positively identify every detected 
feature within the ±600 km s~^ range at the position of O VI. Most such features not associated 
with the targeted galaxy are positively identified HI or metal lines from unaffiUated absorbers 
at other redshifts along the sightline. We generally confirm the presence of an O VI absorber 
by detecting the doublet in the expected 2:1 absorption ratio, by excluding alternative identifi- 
cations for the line (such as another ionization species from another cataloged absorption-line 
system along the sightline), and by requiring that the O VI is aligned to within reasonable limits 
with other species in its absorption system. In the few cases where the data do not cover both 
lines of the doublet because of detector edges or contaminating Galactic absorption, the doublet 
ratio is not available but the latter two criteria must be satisfied to count as a detection. 

Most contaminating absorption that directly impacts the detected O VI is easily accounted 
for in profile fits as nuisance absorption because it affects one line of the doublet but not both. 
After all these steps to identify and/or account for contaminating absorption, we conclude that 
any residual contamination to the detected O VI is of negligible impact to our statistical results. 
We confirm non-detections or "misses" with at least one line of the doublet. For absorbers with 
both line of the doublet detected and uncontaminated by unrelated absorption, we have found 
that the optical depths of the two profiles agree to within errors over the full range of velocity, 
which indicates that partial covering of the background source and Une saturation are not im- 
portant effects on our measurements. 

S4 O VI Column Densities and Kinematics 

We obtain column density and kinematic estimates using both direct integration and profile- 
fitting techniques, which complement one another and reduce total systematic error. Each line 
of each O VI absorber has its total column density A^ovi and the corresponding errors obtained 
by direct integration (55) over velocity limits that were visually determined for each system and 
applied to both lines of the doublet. For non-detections, 3a upper limits to the line strength are 
obtained from estimates of the RMS noise in the continuum where the line is not detected, typi- 
cally over the range ±50 km s~^ around the corresponding Lyman series lines of the absorption 
system, or around the galaxy systemic redshift, in cases where no Lya absorption is detected. 

Profile fitting improves on direct line integration of column densities using information 
about the line shapes and placement to constrain the fit. It also helps assess line saturation by 
taking into account the true Une-spread function (LSF) of the instrument. While fitting can give 
good estimates of column density, intrinsic linewidth, and component structure, the derived fits 
can be degenerate and/or non-unique in cases with closely blended components. Nevertheless 
the systematic effects that complicate profile fitting are different from those that affect direct 
integration, so we use both techniques to assess line saturation and minimize total systematic 
errors. 

Our profile-fitting procedure goes as follows. We extract from each reduced COS spectrum 
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two short slices of ±600 km s~^ around each target systemic redshift, at the position of the two 
O VI lines. The line profiles are fitted using Voigt profiles with parameters A^, v, and b obtained 
by minimization. Most well-detected O VI systems have the fit parameters tied together 
for both lines of the doublet; this tying is also done separately for each velocity component 
where there is more than one. The lines are not tied together where one is unavailable or 
too contaminated to lead to a good fit. Nuisance and/or contaminating lines are fitted where 
necessary (usually with arbitrary atomic data). The O VI absorbers generally require 1 — 4 
velocity components to generate a good fit, with roughly two-thirds of the sample showing only 
one component at COS resolution, about one-fifth showing two components, and the remaining 
systems three or four components. The fitting uses the actual COS LSF (56) based on the fitted 
line's observed wavelength and convolves the model spectrum with this LSF before comparing 
the model to the data. The result of this fitting process is a set of fitted parameters (A^, h, v) for 
each component in each system, and their errors derived from the parameter covariance matrix 
around the best-fit point. The statistical errors on iVovi are generally < 0.1 dex. We find that 
the saturation corrections implied by the difference between direct integration and the profile- 
fitted column densities are small, with the fitted values generally being larger by < 0.1 — 0.2 
dex (which also lends confidence to the errors estimated from covariance of the fit parameters). 
We adopt the fitted values in the main analysis. If the saturation correction is underestimated at 
COS resolution then the A^ovi generally increase and our main conclusions get stronger. 

The strong O VI we detect is in close velocity association with the galaxies. Figure 2 of the 
main text shows the component velocity centroids and the full range of the O VI absorption with 
respect to the galaxy systemic redshift. The errors on the component velocities are generally 
±30 km s^^ and arise from a combination of the 25 km s^^ in the galaxy systemic redshift 
and 10 — 15 km s~^ errors in the COS wavelength solution; multiple components in the same 
system would shift together owing to errors in either one of these quantities. The distribution of 
the detected component velocities has a mean of 18 km s"^ and a standard deviation of 108 km 
s~^, and so the strong O VI we detect is closely associated with the sample galaxies in velocity 
space. This is perhaps not unexpected, given that we have deliberately probed close to galaxies, 
but there is nothing in the close physical separation from the sightlines that requires the O VI to 
favor the galaxy systemic velocity, as we find, over the broader range over which we conducted 
our search. 

To compare these velocities with the escape velocities of the halos in which they reside, we 
convert the measured stellar masses to total dark-matter halo masses using the relation of {58). 
We then assume a spherically symmetric NFW profile (c = 15) and calculate the escape velocity 
as a function of halo mass at i? = 50, 100, and 150 kpc (Figure 2 of the main text). The absorp- 
tion does not tend to significantly exceed the estimated escape velocities of the galaxies. This 
is shown by the cyan range bars on the data points in the lower panel, which mark the velocity 
range over which each absorber was integrated (e.g. full width at zero optical depth). Some 
absorbers have their full ranges touch or slightly exceed tiesc, but this is generally attributable to 
line wings. 

Since we observe the detected material projected on the sky, we measure its velocity rel- 
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Figure SI: O VI component velocities relative to local escape velocity. The dotted line marks 
a ratio of unity, corresponding to the local escape velocity, and the dot-dashed line marks the 
relative galaxy/absorber velocity divided by y/S times the local escape velocity. 

ative to the galaxy projected onto the line of sight instead of its true space motion. Thus the 
measured relative velocities are strictly lower limits to the relative velocity between galaxy and 
absorbing cloud. To assess the importance of this effect, we show in Figure SI the relative 
galaxy/absorber velocity data shown in Figure 2 of the main text, now divided by its local es- 
cape velocity (for its R). Thus any material with a line-of- sight velocity exactly equal to fesc 
will fall on the dotted lines at ±1. To account for unconstrained projection effects, we plot 
also the relative galaxy /absorber velocity divided by v^i^csc^ where the prefactor comes from 
the mean correction between three- space velocity and the line-of- sight velocity for a random 
velocity field. Even with this correction, most of the detected absorption falls within the es- 
cape velocity. We therefore conclude that we have no strong evidence that most of the detected 
high-ionization gas is escaping or has escaped its host galaxy halo. 

There remains the chance that these galaxies have associated O VI absorption at higher rel- 
ative velocities that remains undetected because it falls below our S/N limits. Roughly, we can 
estimate that this absorption would have a factor of ~ 10 — 30 less column density than the 
detected absorption, since the detections typically have column densities that are 10 — 30 x the 
detection limits. 
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S5 O VI and Trends with Star Formation 



We find a sharp difference in tlie detection rate and typical strength of O VI absorption 
surrounding galaxies that are star-forming and galaxies that are passively evolving. This trend 
is clearly evident in Figures 2 and 3 of the main text, where we use specific star formation 
rate as the key variable, because it cleanly divides those galaxies where we have detected star 
formation from those where we have not, at sSFR ^ 10^^ yr^. For the 30 galaxies in which 
we have detected star formation (those with sSFR > 10^^^ yr^^), there are 27 detections of O 
VI (binomial probability P = 0.90 ± 0.07 at 95% confidence), while for the 12 galaxies with 
lower sSFR there are 4 O VI detections (P = 0.33 ± 0.27 at 95% confidence). 

We can also test the statistical significance of the apparent correlation of O VI with sSFR. 
We use a generalized Kendall's tau statistic that can handle doubly censored data (57). First, 
we establish that the correlation of O VI with star formation (sSFR) is strong, when the passive 
galaxies are included. For all 42 sample galaxies, we obtain r = 0.617 (significance P — 
1.32 X 10~^) and so can reject at > 99.99% confidence the null hypothesis that there is no 
correlation between A^ovi and sSFR. This no-correlation null hypothesis is rejected at > 98% 
in the inner two annular bins of R (Figure 1 in the main text), and at 93% over R — 100 — 160 
kpc. Thus there is strong evidence that the correlation with the presence or absence of star 
formation holds at all radii in the detected halos. 

The well-known correlation of galaxy star formation with stellar mass is potentially a con- 
fusing influence on the O VT, one which hinders our ability to conclusively identify star for- 
mation as the key influence on halo O VI. However, even though the size of our sample is 
relatively small, we have evidence that star formation itself, and not simply stellar mass, has a 
direct relationship with the observed O VI. To assess this case we consider the star-forming vs. 
passive comparison in the range of stellar mass where we have both types of galaxies. We take 
subsamples with log > 10.5, where there are 11 star-forming galaxies (8 detections) and 12 
passive galaxies (4 detections). 

We first consider how the quantity of observed O VI differs between star-forming and pas- 
sive galaxies in the same mass range. This comparison appears in Figure S2, where we repeat 
the Nqyi versus sSFR plot shown in the upper panel of Figure 3 in the main text but this time 
including only those galaxies with logM* > 10.5. Here we see that the strong difference in 
hit rates between star-forming and passive galaxies is attributable to the very different column 
density distributions that they follow. The detections of O VI surrounding passive galaxies are 
generally weaker than the detections in the star-forming sample, while the non-detections sur- 
rounding passive galaxies imply weaker absorption than almost all the detections in either cat- 
egory. If we compare the star-forming and passive galaxies in terms of A^ovi using a two-sided 
Kolmogorov-Smirnov test, we find that we can reject at > 99% confidence the null hypothe- 
sis that the two column density subsamples are drawn from the same parent distribution (KS 
statistic D = 0.644, probability 0.00866). 

As a consequence of the lower typical O VI column density surrounding passive galaxies, 
their hit rate is lower than for the star-forming galaxies. We can use the Wilson score interval to 
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estimate the distribution of the underlying binomial probability of detection, given an observed 
sample size and hit rate. For the star- forming galaxies, the underlying binomial probability of 
a detection lies in the range 0.36 — 0.93 at 99% confidence with a hit rate of 0.73. For passive 
galaxies, the 99% confidence interval is 0.36 — 0.68 with a hit rate of 0.33. Thus, the hit rate for 
each category lies just outside the 99% confidence interval on the binomial probability given 
by the other, showing that the two rates are different at ~ 2.6a. Because of the lower column 
densities and hit rates for passive galaxies, even when we control for stellar mass, we conclude 
that the dichotomy of O VI is related to star formation, not solely to mass. We still lack the 
statistics necessary to identify what fraction of the variation in O VI between the blue cloud and 
the red sequence is related to star formation and what fraction owes to stellar mass. 
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Figure S2: O VI vs. sSFR for only those galaxies with stellar masses logM^, > 10.5. This 
figure includes all 12 passive galaxies and 1 1 of the star-forming galaxies from the main sample. 
The column density distributions of O VI are clearly different even over the stellar mass range 
that star-forming and passive galaxies have in common. 

From inspection of Figure 3 in the main text and Figures S2 and S3 here, it is evident that 
much of the statistical power of the full sample resides in the red galaxies with no detected star 
formation and little O VI. We can also test whether the O VI is influenced by the quantity of 
star formation, or only its presence or absence. To do this we omit the passive galaxies from the 
sample and use only the 30 galaxies for which star formation was detected. 

Removing the SFR upper limits weakens these correlations significantly. For the three 50 
kpc annular bins the no-correlation null hypothesis can only be rejected at ~ 75 — 80 % confi- 
dence (Spearman's p ~ 0.37 — 0.39, P ~ 0.21 — 0.26). However, for the whole dataset over 
R = — 160 kpc, we still see an indication (p = 0.352, P = 94% confidence with = 30) of 
a correlation of O VI with sSFR. Though these correlations are less statistically significant than 
those for the whole sample including passive galaxies, there is an indication that the quantity of 
O VI is at least weakly dependent on the quantity of nearby star formation and not simply to its 
presence or absence. 



22 



E 
o 



> 
o 



O 



15.5 



15.0 



14.5 



13.0 



□ Star-forming galaxies 
O ^ Passive galaxies 



°o CD in da 



4a ^ '"^ 



14.0 h ♦ ^ 



o 



< R < 160 

I i_ 



0.01 0.10 1.00 10.00 

SFR yr-^] 



Figure S3: O VI vs. SFR over all R in the sample out to 160 kpc. 



Previous characterizations of gas in galaxy halos at low redshift have not found a strong 
correlation between halo gas and galaxy bimodality. The behavior of O VI in halos contrasts 
with that of the cold gas clouds traced by singly-ionized magnesium, Mg II (30), which is rarely 
detected at i? > 100 kpc, has a steeper radial dependence within that radius, and has a positive 
correlation with galaxy luminosity that is opposite to our findings for O VI. Our findings on 
O VI are in better accord with the 140 kpc size of the triply-ionized carbon (C IV) envelopes 
surrounding galaxies as reported by (59), but that study did not detect a correlation with star- 
formation (as proxied by morphological type). Thus the O VI, C IV, and Mg II may trace 
distinct phases of the circumgalactic medium with possibly different origins. Studies of star- 
forming "Lyman-break" galaxies at 5; ~ 2 — 3 have found a connection between star-formation 
and outflowing material traced by Lya and low and intermediate ions (6, 7), and total masses of 
cold gas comparable to that in the galaxies themselves. 

S6 Physical Conditions, Distribution, and Mass of the O VI- 
bearing Halo Gas 

The narrow range of physical conditions in diffuse gas that contains significant O VI allows 
us to obtain robust lower limits to the mass of oxygen in the detected halo gas. O VI is not 
the dominant ion of oxygen in astrophysical conditions, so the ionization correction applied to 
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derive the total oxygen density (or mass) from the O VI density (or mass) is always substantial 
(Figure 4 of the main text). To estimate the correction, we have considered the ionization 
state of O VI over a wide range of temperature as computed by the CLOUDY photoionization 
code (60) assuming ionization equilibrium and including both coUisional ionization and photo- 
ionization with the z — 0.2 radiation background from (61). Under plausible CGM physical 
conditions, the ionization timescales are usually 10^ — 10^ years or less, compared to halo 
dynamical timescales of ~ 10^ years, so it seems unlikely that a large fraction of the CGM gas 
could be far from ionization equilibrium. 

To calculate the mass implied by our detection of strong O VI surrounding star-forming 
galaxies, we can exploit the fact that column densities are not sensitive to the exact distribution 
of densities along the line of sight. We can therefore apply basic geometry to deduce the total 
mass of O VI in a volume projecting a circular region of R = 150 kpc on the sky with a typical 
column density log A^ovi = 14.5. The total mass of O VI in solar masses is: 

Movi = rfA^ovi(16m,j). (3) 

Scaling to typical values for our sample and including the O VI ionization correction, we 
obtain the total oxygen mass: 

which has the same functional form and scalings as the relation in the main text. As in the main 
text we also correct the coefficient from 1.4 down to 1.2 according to the measured covering 

fraction around star-forming galaxies in three radial bins: 11/11 detections with i? < 50 kpc, 
9 of 10 with 50 < R < 100, and 7/9 with 100 < < 150. Taking the mean column density 
separately in each of these three bins gives log(iVovi) = 14.7, 14.6, and 14.5 respectively. 
Adopting the (A^ovi) and fhu for the individual bins gives a mass coefficient 1.4 x lO^M©, 
but to impose a firm lower limit for the main text we have adopted the lower envelope of the 
star-forming detections at log(A^ovi) — 14.5 to avoid being skewed by the few high values at 
small R. 

The interstellar gas masses cited in the text were derived from the individual stellar masses 
by applying the mean gas fraction versus stellar mass relation as tabulated by (27). The typical 
interstellar oxygen masses in the galaxies were then obtained by applying the metalUcities, 
Z = 0.5 — 1.2^0, obtained from the measurements of nebular emission lines. These values 
define the green band in Figure 4 of the main text. 
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Figure S4: The ratio between our inferred minimum CGM oxygen mass (Mq = 1.2 x 
IO^Mq, see Equation 2) and the typical ISM oxygen mass Mjgj^, as a function of galaxy stellar 
mass, and with the conservative O VI ionization correction /ovi = 0.2. The blue curve and band 
correspond to the mean observed MZR and its systematic uncertainties. The CGM oxygen mass 
Mo range from 70% of M^-^ at the lower end of the mass range in our sample, IO^'^Mq, down 
to 10% at upper mass limit of IO^^Mq. Because the CGM oxygen mass increases in inverse 
proportion to fovi, the blue band expresses a lower limit to the CGM-to-ISM ratio. The heavy 
dashed line shows the ratio between the inferred minimum CGM mass and the oxygen mass 
produced by these galaxies over the last 1 .5 Gyr. 

We have also calculated Mf^^.j for typical star-forming galaxies, using a mean gas fraction 
and mean mass-metallicity relation for SDSS galaxies following the analysis of (27). These 
results are shown in Figures 4 and S4. Here we have derived the gas fractions of galaxies 
using the total budget of cold gas (atomic plus molecular). The mean ISM oxygen mass as a 
function of galaxy stellar mass is then derived and compared with the minimum halo oxygen 
mass, 1.2 X 10'^ Mq x (0.2//ovi)- The widths of these hashed regions express the systematic 
uncertainty in the mass-metallicity relation. For the mean mass-metallicity relation adopting 
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the total cold gas budget, we find that Mq /M^y^ = 0.7 at log M* = 9.5 and 0.1 at log = 11 
(Figure S4). This ratio increases as /ovi declines from its maximal value of 0.2. Notably, the 
ratio between CGM and ISM oxygen masses increase as galaxy stellar mass declines, if /ovi is 
held fixed. 

We can also consider the ratio between the measured CGM oxygen mass Mq and the oxygen 
produced by star formation over some timescale. In this comparison we adopt the mean SDSS 
SFRs as a function of stellar mass, the oxygen yield of 0.014 Mq per Mq of star formation, 
and a timescale of 1.5 Gyr, the travel time to 150 kpc of a cloud moving at 100 km s~^. The 
minimum ratio between the detected CGM oxygen mass and this "produced" mass is given by 
the heavy dashed line in Figure S4. As stated in the main text, the detected CGM mass is a 
notable fraction of the total amount of oxygen produced by these galaxies over a significant 
timescale. 

Because we measure column densities, many possible arrangements of gas along the line 
of sight are possible. However, we can estimate the total line-of-sight extent of the gas using 
simple scalings. For a solar (Z = Zq) oxygen abundance (26), the O VI column density is 

where L^pc is the path length in kpc and p is the mean cosmic density at z = 0.2 (particle density 
of uh = 3.26 X 10^^ cm^^). Thus halo clouds with density uh = 10^^ cm^^ and /ovi = 0.2 
can give logiVovi = 14.5 if they extend 10 kpc along the line of sight, but the required path 
length rises to 200 kpc if the ionization fraction of O VI is only 1% rather than its peak value 
of 20% (Figure S5). Except at very low densities, where the cloud sizes become too large to 
fit within galaxy halos and still give the column densities we observe for star-forming galaxies, 
the ionization fraction of O VI is never more than /ovi = 0.2. 

The cited oxygen masses for the Mg II halos of galaxies were obtained by taking the total 
baryonic mass estimate for that medium from (30), correcting down by 0.1 solar metallicity 
as argued in that paper, and then applying the solar relative abundances of magnesium and 
oxygen. The cited mass estimate for Galactic HVCs is taken from the estimate by (28) for the 
more conservative scenario that locates all the HVCs within 60 kpc of the Milky Way disk. 

To calculate the fraction of cosmic oxygen represented by the ionized halos of star forming 
galaxies, we apply the generic R — 150 kpc halo with A^ovi = 10^^-^ cm~^ to all star-forming 
galaxies with M* > 10^-^ Mq, the range covered by our survey. Integrating over the repre- 
sentative population (specified by the K-selected late-type luminosity function of (50), with 
a = -0.94), we obtain Vl^^^f = 8.2Q/h x 10"^ x (l//ovi)- The cosmic density of oxy- 
gen is rio = yVt^ = 0.014 x 0.00197//i = 2.76/h x 10"^ The fraction of cosmic metals 
in star-forming halos like those in our sample is then f2ov?/^o — 0.015 x (0.2//ovi), as 
given in the text. If we correct this cosmic mass budget to the total gas mass implied in these 
halos, scaling to the Milky Way HVC abundances, we obtain Q^^^"^ = 177{ZQ/Z)n^^f = 
7.31/h X 10~^ X (O.IZq/Z) X (0.2//ovi)- Dividing by the cosmic baryon fraction (52), we 
obtain nl""^"' = O.O23(Z0/Z) x (0.2//ovi). Thus the ionized gas in the extended halos of our 
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star forming sample represent at least 0.2% of cosmic baryons if they have /ovi = 0.2 and solar 
metallicity, 2% if conditions are not optimized for O VI (/ovi ~ 02), or if the metallicity is 
O.IZq, or up to 20% if they have both low metallicity and low /ovi- The contribution of these 
halos to the baryon budget is similar to or exceeds that from atomic and molecular phases of 
the ISM (62), which also make 1 — 2% contributions, but only if conditions are not optimal for 
O VI and/or the gas metallicity is low. 

Figure S5 suggests that it is just possible to explain our observed column densities with solar 
metallicity, photo-ionized gas at T < 10^ K, provided that this gas fills enough of the halo to 
accumulate path lengths > 100 kpc. However, for Z = O.IZq the required path lengths become 
much larger than the virial radii of the galaxy halos, expected to be 200 — 350 kpc for galaxies 
in our observed stellar mass range. Furthermore, if halo gas traces the dark matter, then typical 
overdensities at radii ~ 100 kpc are several hundred or more, in which case photo-ionization 
is ineffective. These arguments favor coUisionally ionized gas, which requires finely tuned 
temperatures to achieve high /ovi (see Figure 4 and discussion in the main text). There are two 
significant caveats to the argument for colhsional ionization. First, the gas responsible for the 
OVI absorption could have lower overdensity than the halo dark matter, with (p/p)gas ~ 50 — 
100 at i? ~ 100 kpc. Large path lengths like those shown in Figure S5 would still be required, so 
this solution is only viable if the metallicity is close to solar. Second, the UV background at the 
10 Rydberg ionization threshold of O VI could be a factor of several higher than the value we 
have adopted (61), in which case photo-ionization becomes important at an overdensity higher 
by the same factor. To give a concrete example, if the UV background intensity is a factor of four 
higher, then /ovi ~ 0.05 at T = 10^-^ K and p/p = 400 (see the p/p = 100 curve in Figure 4), 
and the required path length for Nqyi = 10^^-^ cm^^ drops from 100{Zq/Z) kpc (Figure S5, 
red curve) to 25(Zq/Z) kpc. With our current data, we cannot determine the metallicity of the 
absorbing gas or decide between collisionally ionized and photo-ionized scenarios (or a mix 
of the two). In any case, however, achieving our high covering factors and observed column 
densities requires that favorable combinations of physical conditions arise in the typical halos 
of star-forming galaxies. 
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Figure S5: Path length required to obtain Nqyi = 10^^-^ cmr'^ for solar metallicity gas 
with the indicated overdensity and temperature, using the /ovi values shown in Figure 4 of the 
main text. These pathlengths scale in inverse proportion to metallicity according to Equation 
(3), and so are 10 x longer for Z = O.IZq. The overdensity and ionizing background are 
computed at z = 0.2, using the UV background model of (61). Gas with p/p > 1000 is 
close to coUisional ionization equilibrium, so /ovi depends only on temperature and required 
path lengths are inversely proportional to density. At overdensity p/p < 100, photo-ionization 
becomes important for T < 10^^ K, so /ovi depends on the gas density and the scaling of 
required path length with density is no longer simple. 
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